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Abstract 

The dynamics of vortices in a 2D Heisenberg antiferromagnet with an easy-plane 
anisotropy is studied numericaUy within the discrete spin model as well as analytically 
within a continuum approximation based on a suitable extension of the relativistic 
nonlinear cr model. We find that two like vortices scatter at 90° during a head- 
on collision, whereas a vortex-antivortex pair is annihilated into spinwave radiation 
emitted mainly at 90°. When a uniform bias field is applied, vortex dynamics is 
affected rather profoundly and acquires the characteristic features of the Hall effect 
of electrodynamics or the Magnus effect of fiuid dynamics. In particular, a single vortex 
is always spontaneously pinned, two like vortices form a rotating bound state, and a 
vortex-antivortex pair undergoes Kelvin motion. Finally, in the presence of a bias field, 
vortices are shown to be the prominent topological excitations even for an isotropic 
antiferromagnet . 



I. INTRODUCTION 



Topological magnetic solitons have been studied extensively in the case of 
ferromagnets (FM) and weak ferromagnets (WFM). In both cases a nonvanishing 
magnetization develops in the ground state, albeit by a different physical mechanism, 
which allows a detailed experimental investigation by standard techniques [1,2]. In 
contrast, direct experimental evidence for pure antiferromagnetic (AFM) solitons is 
limited. Nevertheless, theoretical arguments suggest that static AFM solitons should 
exist for essentially the same reason as in ordinary ferromagnets, even though their 
dynamics is significantly diff'erent. 

The dynamics is now governed by suitable extensions of the relativistic nonlinear 
(7 model [2] instead of the Landau-Lifshitz equation [1]. Therefore dynamical concepts 
familiar from the theory of FM domain walls and bubbles need to be reanalyzed within 
the relativistic theory. It is the general purpose of this paper to pursue a study of the 
dynamics of two-dimensional (2D) AFM solitons that emphasizes the influence of the 
underlying topological structure. Our starting point is some recent work on AFM and 
WFM domain walls [3] which revealed that some important issues in the derivation of 
the associated nonlinear a model had been mistreated in earlier treatments. We were 
thus sufficiently motivated to extend the analysis of Refs. [3] to the case of layered or 
2D antiferromagnets whose significance has increased in recent years in connection with 
high-Tc superconductivity. 

The present work proceeds by a combination of numerical and analytical methods. 
Numerical calculations were performed within the standard discrete spin model. 
However a transparent interpretation of the numerical results could not be achieved 
without the aid of a continuum approximation at the heart of which lies the relativistic 
nonlinear cr model. A complete account of the continuum model is given in Sec. II 
where we discuss, in particular, the relevance of certain parity-breaking contributions 
that are implicit in spin models involving antiferromagnetic interactions. Numerical 
and analytical results are then combined in Sec. Ill to provide a complete description 
of static AFM vortices. 

Subsequent sections are devoted to a detailed study of some special dynamical 
features due to the vmdcrlying topology. For instance, head-on collisions of vortices are 
examined in Sec. IV and shown to exhibit a characteristic 90° scattering pattern familiar 
from similar studies of relativistic skyrmions and monopoles [4]. The most surprising 
element is described in Sec. V where it is shown that an applied uniform magnetic field 
affects vortex dynamics rather profoundly. A direct link between topology and dynamics 
is established by means of the conservation laws of linear and angular momentum 
expressed as moments of a suitable topological vorticity, in analogy with related 
work for magnetic bubbles in ordinary ferromagnets [5, 6]. Possible phenomenological 
implications of the derived dynamical picture are further discussed in the concluding 
Sec. VI. Finally, an Appendix presents some useful virial relations that generalize the 
well-known scaling theorem of Derrick [7] . 
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II. THE NONLINEAR (j MODEL 



We shall study a spin system described by the Hamiltonian 

w = ■ (S^+i,. + S^,,+i) + IgiSljf], (2.1) 

where the summation extends over all sites of a square lattice, S^^- is the spin vector 
at site and Sfj its third component. Both the exchange constant J and the 

anisotropy constant g are taken to be positive and hence (2.1) describes a Heisenberg 
antiferromagnet with an easy-plane single ion anisotropy. 

The spin variables are treated as classical vectors with constant magnitude s and 
satisfy the equation of motion 

where the effective field F is determined from the general relation 

dW 

^'■^ = -as-- P-^' 

or, more explicitly, by 

Fij = -J{Si+i,j + Sj-ij + Sij+i + Sij-i) - gSfjB, (2.4) 

where e = (0, 0, 1) is a unit vector along the third direction. The exchange contribution 
in the effective field (2.4) contains fewer than four terms on the perimeter of a finite 
lattice, the precise number of such terms being equal to the number of nearest neighbors. 
All dynamical simulations presented in this paper will be based on the above relatively 
simple set of discrete equations adapted to an open finite lattice. 

However, in order to pursue an efficient study of the dynamics, we also sonsider a 
continuum approximation which is possible at weak anisotropy, 

(2.5) 

where soliton structures extend over a large number of sites given roughly by 
1/s. The continuum limit is not completely straightforward because of the implicit 
antiferromagnetic discontinuity as one moves from site to site. It is then important to 
first identify dynamical variables that may possess a smooth limit as e ^ 0. 

In one dimension continuity is achieved by a simple dimerization process [3]. A 
square lattice may also be thought of as a collection of dimers, as illustrated in Fig. 1 
for a finite lattice cut along the diagonals. The sites of the original square lattice are 
depicted by two sets of (solid and open) circles which form two intertwining sublattices. 
This designation is used merely to indicate that the two spins of any given dimer 
belong to different sublattices. In particular, the ground (Neel) state is such that the 
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two spins point in opposite directions but are uniform on each sublattice. Because of 
the easy-plane anisotropy the Neel state is polarized along any direction in the plane 
perpendicular to the third axis. This azimuthal degeneracy of the ground state will play 
an important role in the following. 

A generic dimer AB (see Fig. 1) is labeled by a pair of indices (a, /3) numbered 
consecutively {a, P = 1,2, ... , N). In our explicit illustrations N is taken to be even, 
but this is only a minor technical assumption. Let us denote by Aa^p and Bq,,/3 the two 
spins of dimer AB. Eq. (2.2) is then written as a system of two coupled equations, 

where the effective fields F and G are given by 

Fa,/3 = -</(Ba,/3 + B^-i^/j + Ba,/3-l + Bq,_i,^_i) - gA^^ j^^, 



(2.7) 



for a generic point inside the lattice. The exchange contributions in Eq. (2.7) contain 
only two terms for each point on the perimeter of the finite lattice of Fig. 1 . 

As a first step in the derivation of a continuum approximation we introduce the 
discrete set of variables 

r] = V2e{a-ao), ^ = V2e{p-Po), (2.8) 

which become continuous in the limit £ — > and provide a measure of distances along 
the diagonals of the square lattice. Actual distances are given by arj/s and a^/e where 
a is the physical distance between two neighboring magnetic ions. The lattice constant 
a will not appear in our theoretical development except when various quantities of 
interest will have to be translated in physical units. Finally the convenient choice 
ao = {N + l)/2 = /?o in Eq. (2.8) sets the origin of the coordinate system at the center 
of the lattice of Fig. 1. 

The main assumption supported by numerical calculations is that the spin variables 
Aa^l3 and Bq,^^ approach smooth continuum limits A = A(?7, ^) and B = B(77, ^) at weak 
anisotropy (e — > 0). Then we make the substitutions A^t^p — > A and Bq,^^ ^ B in Eqs. 
(2.6)-(2.7) together with the Taylor expansions 

^a±l,/3 ^ A ± SArj + ^S'^Ar^rj, 

A,,^±i ^A±5A^ + Is^A^^, (2.9) 

A ± S{Arj + A^) + -S'^iAr^n + A^^ + 2A^^), 

and similar expansions for the field B. Here subscripts denote differentiation with 
respect to the indicated arguments and 

S = V2e (2.10) 
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is used as a temporary notational abbreviation. Eqs. (2.6) are then approximated by 

^ = AxF, ^ = BxG, (2.11) 

where 

F = - J[4B - 2S(Br, + B^) + S'^iBrjrj + B^^ + B^^)] - gAse, 

(2.12) 

G = - J[4A + 25(A^ + Ae) + ^'(A^^ + A^^ + A^^)] - ^^36, 

This system of equations is not yet fuUy consistent because it appears to mix different 
powers of S (or e). 

In order to obtain a consistent continuum model we proceed as in the ID case 
studied in Ref. [3]. First, we introduce the hnear combination of fields 

m=-^(A + B), n=l(A-B), (2.13) 

which satisfy the constraints 

m-n = 0, m^ + n^ = l. (2.14) 
Second, we define a dimensionless time variable 

T = 26sJt = 2V2esJt. (2.15) 
An equivalent form of Eqs. (2.11) then reads 

1 

5— = - 5[{m. X n)^ + (m x n)^] + -5^[n x (n^^ + n^^ + n^^) 

- m X (m^^ + m^^ + m^^)] - [777,3(111 x e) + 773(11 x e)]. 



(2.16) 



dn 

S— = 4(m x n) + 5[m x (m^ + m^) — n x (n^ + n^)] 

— -^5^ [7773(11 X e) + 713(111 X e)]. 

A simple inspection of the above equations suggests that consistency is obtained if m 
is of order S. Then the leading approximation of the second equation is 

S— = 4(m X n) - 6[n x (n^ + n^)] (2.17) 

and the constraints of Eq. (2.14) reduce to 

m-n = 0, 11^ = 1, (2.18) 
5 



to within terms of order 5"^. Therefore taking the cross product of both sides of Eq. 
(2.17) with n and using the constraints (2.18) yields 



m = ^[-K + n^) + (n X h)], (2.19) 

where we have restored the small parameter e from Eq. (2.10). Finally Eq. (2.19) is 
inserted in the first of Eqs. (2.16) to give in the limit £ — > the differential equation 

nxf=0, f=h-An + n3e, (2.20) 

where A is the 2D Laplacian 

" +1^- (2-21) 



It is understood that terms of order £^ in Eqs. (2.18), (2.20) and order in Eq. (2.19) 
have been neglected. 

Therefore the continuum approximation is governed mainly by Eq. (2.20) which is 
a simple extension of the relativistic nonlinear a model to include a single ion anisotropy. 
The corresponding velocity of light is equal to unity thanks to our choice of rationalized 
space and time variables. In our conventions, the spin magnitude s is a simple multiple of 
the Planck constant (e.g., s = ^h) and thus carries dimension of action, sJ of frequency 
and J of energy. The ratio e of Eq. (2.5) is dimensionless, so is the time variable r of 
Eq. (2.15). Also recalling that actual distances are given by ar]/e and a^/s, where a is 
the lattice constant, we conclude that velocity is measured in units of 

c = 2V2asJ, (2.22) 

which coincides with the phase velocity of pure AFM magnons on a square lattice in 
the long- wavelength limit. More generally, the magnon velocity is given by 

c = 2asJ^, (2.23) 

where z is the lattice coordination number. Applied for a square lattice {z = 4) Eq. 
(2.23) reduces to Eq. (2.22), whereas for a chain {z — 2) it yields the limiting velocity 
c = 2as J of Refs. [3]. 

Needless to say, a complete description of the original spin model requires knowledge 
of both fields m and n. On the other hand, given a solution n = n(?7,^,r) of 
the nonlinear a model, the field m may be determined from Eq. (2.19) by simple 
differentiations and may thus be viewed as an auxiliary field. Yet this apparently 
straightforward result has been controversial. For instance, the parity-breaking gradient 
terms in Eq. (2.19) were recently derived within the ID model and shown to be important 
for various structural properties of AFM and WFM domain walls [3]. The possible 
occurrence of such terms had been anticipated on symmetry grounds [8] but this 
possibility was overlooked in the literature for a long time [9]. 
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One should stress that symmetry arguments do not predict the precise coefficients 
of the parity-breaking gradient terms in Eq. (2.19) and arc generally susceptible to 
overinterpretation. The danger from this dimerization ambiguity is already present 
in one dimension and was completely analyzed in Refs. [3]. The situation is only 
compounded in two dimensions. For example, had we chosen to work with a horizontal 
instead of the vertical dimerization of Fig. 1, the nonlinear a model of Eq. (2.20) would 
not be affected but the gradient terms in Eq. (2.19) would appear as — instead 
of + n^. This is an indication that the local values of the field m are sensitive to the 
mode of dimerization and cannot be literally interpreted as magnetization. However no 
real mathematical or physical ambiguity appears when the results of a calculation or an 
experiment are consistently interpreted in reference to a specific mode of dimerization. 
But an attempt to measure the "magnetization" m by standard techniques may lead to 
a fuzzy magnetization curve around any nontrivial soliton structure, unless spin values 
can be resolved at every site [3]. 

To press the above picture further we note that dimerization is not the only way 
to achieve continuity on a square antifcrromagnetic lattice. We may also consider the 
tetramerous configuration illustrated for a finite lattice in Fig. 2. Each tetramer is 
again labeled by two indices a and ^ numbered consecutively (a, ^ — 1,2,..., A^) and 
spin values are smooth as one moves from corresponding sites of one tetramer to the 
next. Let us again denote by Aq,^^, Bq,^^, C^^/j and D^,/? the spins on a generic tetramer 
ABCD shown in Fig. 2, in terms of which the original equation of motion (2.2) may be 
written as a system of four coupled equations analogous to Eqs. (2.6). This new system 
may then be used for the derivation of the continuum approximation by a method 
similar to the one explained earlier within the dimerization scheme of Fig. 1. 

We omit the lengthy algebraic details and simply state the final results. We now 
use the Cartesian coordinates 



which are related to the coordinates rj and ^ of Eq. (2.8) by a 45° rotation. Here the 
origin of the coordinate system is again set at the center of the lattice of Fig. 2 by 
choosing the arbitrary constants as ao = [N + 1)/2 = Pq. Next we assume that the four 
spins on a tetramer approach smooth continuum limits when £ — > denoted by A, B, C 
and D which are some functions of the spatial coordinates x and y of Eq. (2.24) and 
the time variable r of Eq. (2.15). A more convenient set of fields is given by 



X — 2e{a — ao), 



(2.24) 



m= -^(A + B + C + D), 
4s 

k=-^-(A + B-C-D), 





1 



(A 



B + C-D), 



B-C + D), 



(2.25) 



and satisfy the constraints 



m • k + n • 1 = 0, 



(2.26) 



m • n + k • 1 = 0, 



m • 1 + n • k = 0. 
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In the strict continuum limit (e — > 0) the constraints reduce to 

m-n = k- n = l- n = 0, n^ = l, (2.27) 
the fields m, k and 1 are expressed in terms of n by 

m=^(nxri), k=-|n„ l=-^ny, (2.28) 

and n itself satisfies the differential equation 

nxf=0, f=h-An + rise, (2.29) 
where A is the 2D Laplacian 

It is again understood that Eqs. (2.27)-(2.29) are accurate to within terms of order 

Although the field n defined in Eq. (2.25) is not directly related to the field n of 
Eq. (2.13), its continuum dynamics is still governed by the nonlinear a model of Eq. 
(2.29) which is equivalent to Eq. (2.20). However a clear distinction between the two 
formulations emerges at the level of the auxiliary fields. We now need three such fields 
(m, k and 1) which are expressed in terms of n through Eqs. (2.28). These relations 
reenforce our earlier remarks concerning parity- breaking contributions, namely that 
they are sensitive to the specific mode of taking the continuum limit. For instance, such 
contributions are no longer present in the field m but their effect is accounted for by 
the new auxiliary fields k and 1. Once again we must conclude that the local values of 
the field m cannot be interpreted literally as magnetization. 

The perplexing nature of the auxiliary fields calls for a summary of our main 
strategy. The continuum model is extensively used to motivate the results discussed 
in subsequent sections. But all numerical calculations are based on the original discrete 
equation (2.2). The actual calculations may be performed on an open finite lattice 
either of the type shown in Fig. 1 or that of Fig. 2. The explicit results may then 
be employed to construct the fields m and n in the former case and the fields m, n, k 
and 1 in the latter. In both cases the validity of the respective continuum models can 
be verified explicitly for a sufficiently weak anisotropy {e ^ 1). Nonetheless several 
apparent paradoxes emerged in the course of our investigation which were all resolved 
in favor of the formulation presented in this section. 



III. STATIC VORTICES 

We consider first the problem of finding interesting static solutions within the 
continuum model. In both formulations developed in the previous section, the basic 
issue is to determine the field n from the nonlinear a model and then proceed with the 
calculation of the aiixiliary fields. Time derivatives vanish in the static limit and one 
obtains the reduced equations 

nxf=0, f= -An + 7136. (3.1) 
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It proves convenient to derive the field f from a variational argument, 



SW 1 



6n 



^=\j [(^^n • ^^n) + nl]dxdy, (3.2) 



where W is the energy functional in the static limit (see Sec. V). The repeated index 
is summed over two distinct values jjL — 1 and 2 corresponding to the two spatial 
coordinates x and y (or rj and ^). It is also useful to resolve the constraint = 1 
explicitly using, for example, the spherical parametrization 

ni=sin0cos$, n2 = sin0sin$, ns^cos© (3-3) 

in terms of which 

W=^J [{d^Qd^Q) + sin^ e (d^^d^^) + cos^ Q]dxdy (3.4) 
and Eqs. (3.1) are equivalent to 

or 

AG + [1 - (a^$a^$)] cos e sin e = 0, ^^(sin^ B 5^$) = 0. (3.6) 

It should be noted that the above static equations are formally identical to those 
encountered in easy-plane ferromagnets. However the physical interpretation of the 
field n is now different and the actual construction of the corresponding AFM solitons 
on the discrete lattice is effected by the prescriptions of Sec. II. 

The simplest solution of Eqs. (3.6) is the ground state configuration Q = ^ and 
$ = (po, where (po is an arbitrary constant, for which the energy achieves its absolute 
minimum {W = 0). Such a simple configuration (ni = cos^oi ^2 = sin^o? ''^s = 0) is 
consistent with our earlier remark that the ground state is the usual Neel state with 
spins polarized along any direction in the (12) plane. Nontrivial static solutions are also 
stationary points of the energy functional W and are subject to limitations imposed by 
the familiar scaling theorem of Derrick [7]. Applied for the functional (3.4) Derrick's 
argument yields the virial relation 

cos'^ e dxdy ^0. (3.7) 

Therefore one must conclude that nontrivial static solutions with finite energy do not 
exist. 

The above argument does not exclude interesting solutions with infinite energy 
such as vortices. In fact, the possible existence of vortices is probed by a suitable 
generalization of the Derrick theorem worked out in the Appendix. Eq. (3.7) is then 
replaced by 

cos^ dxdy = tt, (3.^ 
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which not only does not exclude the possibility of vortex solutions but predicts the 
actual value of their anisotropy energy. However the exchange energy is logarithmically 
divergent. 

To obtain an explicit solution we consider the usual cylindrical coordinates x = 
p cos (f), y = p sin and make the following ansatz for a vortex located at the origin: 



Q = 9{p), $ = «(0 + 0o), (3.9) 

where k = ±1 will be referred to as the vortex number and 0o is an arbitrary constant 
reflecting the azimuthal symmetry. Then the second equation in (3.6) is automatically 
satisfled and the first yields 



id_ / de\ r_ 1 

pdp\dp) V P 



The energy (3.4) is accordingly reduced to 

1 



2 







7^ H + cos^ d 



\dpj p 



2 



{2tt pdp) (3.11) 



and its stationary points are solutions of Eq. (3.10). This equation is consistent with 
the asymptotic behavior 

e ^ CP, e ^ f (3.12) 

where c\ and Ci are some constants. The actual solution depicted in Fig. 3 was calculated 
numerically via a relaxation algorithm applied to the restricted energy functional (3.11) 
taking into account the boundary values (3.12). 

Furthermore we note that for every solution d of Eq. (3.10) a new solution is 
obtained by the formal replacement 6 ^ tt — 9. Hence the complete result for the field 
11= (tt-i, ?T'2, 'T's) reads 

ni = sin 6* cos [«;((/) + (/)o)], n2 = sin 6* sin [«;((/> + (/)o)], n3 = i^cos6', (3.13) 

where 9 = 9{p) is taken from Fig. 3, while the vortex number k and the polarity v are 
given by 

K = ±l, u = ±l, (3.14) 

taken in any combination. In other words, both the vortex [k = 1) and the antivortex 
(k = —1) come in two varieties {u = ±1). 

The asymptotic behavior (3.12) makes it evident that the centrifugal (second) 
term in the energy (3.11) diverges logarithmically with the size of the system but the 
remaining terms are finite. In fact, the numerically calculated anisotropy contribution 
was found to be in excellent agreement with virial relation (3.8). The weak (logarithmic) 
divergence of the total energy is not necessarily an obstacle to the production of 
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vortices in an antiferromagnet. In this respect, one should recaU that vortices are easily 
produced in a rotating cylindrical bucket filled with an ordinary fiuid or a superfiuid [10]. 
Furthermore vortex-antivortex pairs have finite energy and should play an important 
role in thermodynamics. 

Keeping with the strategy outlined in the concluding paragraph of Sec. II, we 
must comment on the manner in which Eq. (3.13) furnishes a vortex on the discrete 
lattice. The field n of Eq. (3.13) may be used to calculate the auxiliary field m of 
Eq. (2.19) or the fields m, k and 1 of Eq. (2.28), applied for a static solution (h = 0), 
and thereby determine the original spin configuration on the lattices of Fig. 1 and 
Fig. 2, respectively. As a check of consistency we have used the resulting spin 
configuration as initial condition within a fully-dissipativc algorithm [3] applied to the 
complete energy functional (2.1) on a 200 x 200 lattice. The calculated relaxed state 
accurately reproduced both the profile of Fig. 3 and the corresponding auxiliary fields 
for a reasonably small value of the anisotropy constant {e — 0.1). The spin values on the 
discrete lattice are partially illustrated in Fig. 4 for both a vortex and an antivortex and 
exhibit the characteristic antiferromagnetic discontinuity between the two sublattices. 
Finally we mention that a similar calculation of an AFM vortex can be found in Ref. 
[11] using an Ising-like instead of a single-ion anisotropy. 

Since the vortex solutions described above are present for any finite value of the 
anisotropy, however weak, it is also of interest to examine the extreme limit of vanishing 
anisotropy. The static equations (3.1) then reduce to 

nxAn = 0, = 1 (3.15) 

and are equivalent to the 0(3) nonlinear a model in a 2D Euclidean space. The small 
parameter e hidden in the definition of spatial coordinates is no longer related to an 
anisotropy constant but is an intrinsic scale set by the actual spread of localized solitons. 
The implied arbitrariness of e is reflected by the scale invariance of Eq. (3.15) and is 
already an indication that the corresponding solitons are metastable. Explicit solutions 
of Eq. (3.15) are the well-known Belavin-Polyakov (BP) instantons [12] viewed here as 
static solitons in the 2-1-1 dimensional theory of current interest. They differ from the 
vortex (3.13) mainly in the asymptotic value of the field n which is now given by 

n (0,0,1) (3.16) 

|X|— >oo 

and coincides with the ground state configuration of the isotropic antiferromagnet 
defined up to a global 0(3) rotation. Therefore the BP instantons may be called AFM 
bubbles, by analogy to FM bubbles [1], and are classified by the Pontryagin index or 
winding number Q which is defined as 

If 1 
Q = ^ / qdxdy, q = -£^^.(9^,11 x ^^n) • n, (3.17) 

where e^i, is the 2D antisymmetric tensor, and is integer-valued (Q = 0, ±1, ±2, . . .). 
The topological density q may also be expressed in terms of the spherical variables as 

g = £^^sinea^ea^$ (3.18) 



11 



and plays an important role in the dynamical theory of FM bubbles [5,6]. 

An examination of the asymptotic behavior of Eq. (3.13) suggests that a vortex 
may be viewed roughly as a half bubble. To push this remark further we calculate the 
density q for the vortex configuration (3.13), 

KvdcosO 

<1= 3.19 

P dp 

where 6 = 9{p) must be taken from Fig. 3. Therefore the winding number calculated 
from Eq. (3.17) is found to be 

Q = 1ki/[cos ^(oo) - cos ^(0)] = -^Kiy (3.20) 

and depends on both the vortex numper k and the polarity u. This result confirms 
the vague notion that a vortex is topologically equivalent to a half bubble {Q — ±|). 
However a modified topological charge that is related to the vortex number but not the 
polarity will arise more naturally in the dynamical context of Sec. V. 



IV. HEAD-ON COLLISIONS 

Although static AFM solitons are similar to those encountered in ferromagnets, 
their dynamics is significantly different. For comparison purposes, it is useful to recall at 
this point the two main dynamical features of FM bubbles: (a) An FM bubble cannot be 
found in free translational motion; it is always spontaneously pinned or frozen within the 
ferromagnetic medium, (b) An FM bubble tends to move in a direction perpendicular 
to an applied magnetic field gradient. For example, a single bubble will undergo a 
90° deflection with respect to an externally supplied uniform gradient in the absence 
of dissipation [5, 6]. In the case of two or more interacting bubbles, a gradient arises 
intrinsically and leads to a characteristic relative motion similar to the Hall motion of 
electric charges in a uniform magnetic field or the motion of vortices in a fluid [13]. 

Property (a) is clearly not the case for AFM bubbles or vortices because of the 
Lorentz invariance of the underlying nonlinear a model. Indeed, for any static vortex 
n = n(x) constructed in the preceding section, a vortex propagating freely with an 
arbitrary speed v < 1(= c) in, say, the a;-direction may be obtained by the elementary 
Lorentz transformation 

One may then calculate the corresponding auxiliary flelds from Eq. (2.19) or (2.28), 
taking into account that the dynamical contribution (n x h) is no longer vanishing, and 
subsequently construct the spin conflguration on the discrete lattice to obtain a rigidly 
moving AFM vortex. 

Similarly, there is no reason to believe that property (b) is sustained for AFM 
bubbles or vortices. Yet one should expect that topology will continue to play an 
important role within the relativistic dynamics. A possible manifestation of a topological 
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effect occurs in tlie dynamics of two bubbles in a head-on collision. Such a process was 
studied extensively in the context of the isotropic 0(3) nonlinear a model and shown 
to exhibit a characteristic 90° scattering pattern [14]. This pattern is certainly unusual 
from the point of view of ordinary particle dynamics but strongly reminiscent of the 90° 
deflection of FM bubbles in a field gradient. Of course, the above analogy is superficial 
but indicates that both the Landau-Lifshitz and the relativistic dynamics are influenced 
by the underlying topology. 

The isotropic nonlinear a model is special in several ways. In particular, its scale 
invariance leads to metastable bubbles of arbitrary radius. The radius of each bubble 
changes during collision and never returns to its initial value [14] . A healthier situation 
arises in the presence of anisotropy which sets a definite scale for the soliton size. Then 
individual solitons may be deformed during scattering but will always bounce back to 
their original shape well after collision. The AFM vortices constructed in Sec. Ill are 
examples of 2D solitons with definite scale and their dynamics will be the subject of the 
remainder of this paper. 

The simplest dynamical experiment is to consider the evolution of a pair of two 
idendical vortices which are initially at rest at a relative distance d. Because vortices are 
extended structures, the above initial configuration is not uniquely defined and generally 
depends on the details of the physical process that brings the two vortices to their initial 
positions. However such details may be important for a short transient period but 
are not expected to significantly infiuence the global properties of the ensuing motion, 
especially when the initial distance d is large. Therefore there exists significant freedom 
in the construction of the initial configuration. The simplest choice is to consider the 
product ansatz 

d \ ^ f d 



Vt{x,y) = VL^\^--,yjn^\^ + -,yj, (4.2) 
where Q. is the complex stereographic variable 

^ = r^±±^ ^ tan(e/2) e^* (4.3) 
1 + ns 

for the vortex pair and fii is the corresponding variable for a single vortex. In view of 
Eq. (3.13) we may write 

= 1 e-(^^+^°), (4.4) 
\ 2 / l + v cos d± 

where 




and 6 = 9{p) is the profile of Fig. 3. In our simulations we made the specific choices 
K = 1, = — 1 and 00 = 0. The product ansatz (4.2) is then written as 

n = /(p+)e^^+/(p_)e^^-, /(p±) = cot^ (4.6) 
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and represents a pair of two identical vortices initially at rest, at a relative distance d 
on the X-axis. 

The remaining steps for a complete specification of the initial configuration on the 
discrete lattice of, say. Fig. 2 proceed as follows. The field n is obtained by inverting 
Eq. (4.3), i.e., 

n + n 1 n-U i-oo 

77-1 = ^— , 712 = , n3 = . (4.7) 

The auxiliary fields m, k and 1 are then computed from Eq. (2.28) applied for a static 
configuration (n x h = 0). The original spin variables A,B,C and D on a generic 
tetramer with coordinates {a, (3) are determined from Eq. (2.25) applied for the discrete 
set of points x = 2e{a — ckq) and y — 2e{(3 — (3o) of Eq. (2.24) with a, (3 — 1,2, . . . , N. 
The spin configuration is thus specified at every site of the original lattice. 

Having determined the initial configuration, the ensuing evolution of the vortex 
pair was calculated by a numerical solution of the initial- value problem (2.2)-(2.4) using 
a fourth-order Runge-Kutta algorithm. Typical runs were performed on a 200 x 200 
lattice for a reasonably weak anisotropy e = 0.1 that leads to a respectable grid [—10, 10] 
for the dimensionless position variables x and y. The rigidity of our numerical results 
was frequently checked on larger lattices, up to 400 x 400, and stability was improved 
by reenforcing the constraint S^^- = at every site of the lattice after every Runge- 
Kutta step. This numerical trick was borrowed from related simulations in the isotropic 
nonlinear a model [14]. On this occasion, it should be emphasized that we do not 
directly simulate the dynamics of the a model but rather of the original Heisenberg 
antiferromagnet. A byproduct of this fact was that various theoretical predictions based 
on the continuum approximation of Sec. II were verified in detail. In particular, we have 
been able to illuminate numerous subtle points in connection with the parity-breaking 
contributions in the auxiliary fields. 

The results of the numerical simulation described above revealed no surprises when 
the two vortices are initially at rest; they begin to drift away from each other along the 
X-axis apparently in order to minimize their interaction energy. Therefore the calculated 
behavior is similar to that of two ordinary particles interacting with a repulsive potential. 
Nevertheless this behavior is already significantly difi"erent from the one observed in the 
case of two interacting FM bubbles which would rotate around each other irrespectively 
of whether the potential is repulsive or attractive [13]. 

A more interesting situation arises when the two vortices in Eq. (4.6) are initially 
Lorentz boosted to velocities v and — v and are thus set on a head-on collision course 
on the X-axis. As expected, the two vortices begin to discelerate thanks to their mutual 
repulsion. However the future of the process depends crucially on the magnitude 
of the initial velocity. At low velocities the two vortices approach each other to a 
minimum distance at which they come to rest and then turn around and move off in 
opposite directions. When the initial speed exceeds a certain critical value the vortices 
again discelerate but come sufficiently close to a relative distance where the interaction 
potential seems to have become attractive. The vortices then begin to accelerate toward 
each other until they overlap almost completely. More imporantly, they subsequently 
split and reemerge as two separate vortices moving off along the positive and negative 
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^/-direction, thus undergoing a 90° scattering analogous to the one observed in numerical 
simulations of bubbles in the pure 0(3) model [14]. The main difference here is that 
vortices regain their initial shape except for a flip of their internal phase equal to 
<t>o = ±f . 

The situation described above is illustrated in Fig. 5 which depicts three 
characteristic snapshots of the collision of two like vortices that start from an initial 
relative distance d = 5 with velocities vi = (0.65,0) and V2 = (—0.65,0). We have 
chosen to illustrate only the field n, through its projection on the (12) plane, and hence 
the antiferromagnetic discontinuity present in Fig. 4 is not apparent in Fig. 5. However 
we have actually calculated the spin values on the original lattice and then proceeded 
with the determination of n as well as the auxiliary fields. We were thus able to confirm 
several important details of the formulation presented in Sec. II. 

The picture changes drastically in the case of a vortex-antivortex pair initially at 
rest, at a relative distance d on the x-axis. Such a pair is initially described by the 
product ansatz 

l^ = /(p+)e^'^+/(p_)e-^'^-, (4.8) 

where notation is the same as in Eq. (4.6). The subsequent evolution was again 
studied by a numerical calculation similar to the one described earlier for a pair of 
two like vortices. The vortex and the antivortcx now begin to accelerate toward each 
other apparently because their interaction potential is attractive. Fig. 6 depicts three 
characteristic snapshots of the field n which indicate that the resulting spin configuration 
shows no sign of a nontrivial topological structure, or that the vortex-antivortex pair 
is annihilated. Nevertheless a closer examination of the time evolution of the energy 
density shown in Fig. 7 suggests that the annihilation process is not completely dull. 
Indeed two distinct lumps of energy are emitted along the positive and negative y-axis 
which do not correspond to any well defined soliton structures and eventually dissipate 
into spinwaves. It is interesting that a tendency for a 90° scattering persists even in the 
present case [14]. 

To summarize, our calculation of head-on collisions of AFM vortices confirms what 
appears to be a robust feature of relativistic topological solitons [4, 14]. It is more 
than clear that a simple theoretical explanation of the observed behavior should be 
possible to obtain almost independently of the specific dynamical model. However such 
an explanation has thus far been offered only within the context of some variant or 
other of collective coordinates. For example, vortices in a Higgs model were discussed 
along those lines in Ref. [15] and shown to exhibit 90° scattering, at least in the so-called 
Bogomonly limit. At this point, we do not have a convincing theoretical interpretation 
of our numerical data for head-on collisions. But the theoretical situation improves 
significantly in the case of AFM vortices in a uniform magnetic field studied in the 
following section. 
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V. VORTICES IN A MAGNETIC FIELD 



The effect of an applied uniform magnetic field H is accounted for by including a 
Zeeman term in the Hamiltonian (2.1), 

W^W-gofioJ^C^-^ij)^ (5-1) 

where (70 ~ 2 is the gyromagnetic ratio and /xq = e/2mc is the Bohr magneton divided 
by the Planck constant. AFM vortices then acquire the general dynamical features of 
FM bubbles described in the beginning of Sec. IV. Such a radical change of behavior 
was already speculated in Ref. [16] and is established here by means of unambiguous 
conservation laws that link the dynamics with the underlying topology. The derived 
qualitative picture is then confirmed by direct numerical simulations. 

The classical ground state is now obtained by assigning a spin value A on the 
first sublattice (solid circles) and a value B on the second (open circles). We further 
introduce the unit vectors a = A/s and b = B/s to write for the energy of such a 
configuration 

W/s'jA = 2(a • b) + + bl) - • (a + b), (5.2) 

where s is the dimensionless anisotropy constant of Eq. (2.5) and A is the total number 
of lattice sites assumed to be large. For a field of strength H applied along the third 
direction the minimum of (5.2) is achieved by the canted spin configuration of Fig. 8, 
defined up to an arbitrary azimuthal rotation, where the canting angle is given by 

sm5=—, Hc = , (5.3) 

for field values in the range H < He- Above the critical field a transition takes place 
into a ferromagnetic phase {5 = Actually we shall mostly study the parameter 
regime where both the anisotropy and the applied field are weak, namely 

8s 7 

£<1, H^Hc^ , (5.4) 

which are conditions for the validity of a continuum approximation and are sufficiently 
nonstringent for practical applications. For future reference. Fig. 8 also displays the 
two vectors 

m=-(a + b), n=-(a-b), (5.5) 
which may be expressed in terms of the canting angle as 

m= (0,0, sin 5) ^ ((7o/Uoi7/8sJ)(0,0,l), 

(5.6) 

n = (cos 5, 0,0) fti (1,0,0), 
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where the second steps have been restricted to the parameter regime (5.4). Therefore a 
nonvanishing magnetization m develops along the third direction whereas n is confined 
in the basal plane. 

When the field is turned on, the Neel state is set in a precessional mode that 
eventually relaxes into the canted state of Fig. 8 thanks to some dissipative process or 
other that is always present in a realistic antiferromagnet. Throughout this section 
we shall assume that the field has been turned on sufficiently long to ensure that 
equilibrium has been achieved in the ground state. The argument is carried out in 
five steps described in the following five subsections. 



A. The continuum model 

We now return to the discrete equations (2.6)- (2.7) which we extend according to 
Fa,/3 ^ Fq,^^ + 5fo//oH, Ga,p ^ Ga,/3 + go^ioH- (5.7) 
It is also convenient to introduce the rationalized field [3] 

h = he, h=^^p^. (5.8) 

Eqs. (2.16) are then extended simply by adding a term d{m x h) to the right-hand 

side of the first equation and 5(n x h) to the second, where 5 = \/2e. The remaining 
algebraic details will be omitted here because they are similar to those of Sec. II. 
Hence the continuum limit is now governed by the strong inequalities 

£,e/i-Cl, (5.9) 

the fields m and n continue to satisfy the constraints 

m-n = 0, = 1, (5.10) 

the aiixiliary field m is given by 

m = ^[-("^v + n^) + (n X h) - n X (n X h)] (5.11) 

and the field n satisfies the differential equation 

nxf=0, f= h- An + 2(h X h) + (n-h)h + n3e. (5.12) 

As a check of consistency one may explicitly verify that the ground state values 
(5.6) are compatible with the above equations for parameters that satisfy the strong 
inequalities (5.9). Finally we note that Eq. (5.12) emerges also in connection with the 
tetramerization scheme of Fig. 2, while the associated auxiliary fields are 

m = ^-=[(n X n) — n X (n X h)], 
2v 2 

(5.13) 
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which should be compared to Eqs. (2.28). 

We may thus concentrate on the dynamics of the extended nonhnear cr model 
(5.12) where the effect of the applied field is twofold; it breaks Lorentz invariance and 
renormalizes the single ion anisotropy. An efficient study of the dynamics is carried out 
through a variational principle, namely 

f = -1^, (5.14) 
on 

where A is the action 

^ = y" Ldxdydr, (5.15) 
and L is the corresponding Lagrangian density [2] 

L = ^[h^ - {d^n ■ d^n)] + h • (n X n) - ^[(n • h)2 + n^]. (5.16) 

We further resolve the constraint = 1 explicitly using, for example, the spherical 
parametrization (3.3) to write 

L =^[e^ + sin^ e $2] _ ^[{d^Qd^Q) + sin^ 9 {d^^d^^)] 

+ /i sin^ e $ - ^ (1 + /i^) cos^ e. 
Hence there exist two pairs of canonical fields given by 

■01 = 0, TTi = 

ip2 = ^, 7r2 = sin^0(6 + /i). 
The Hamiltonian is then obtained from 



(5.17) 



(5.18) 



W — J wdxdy, W = TTatpa ~ L, (5.19) 

where the repeated Latin index a is summed over the two distinct values of Eq. (5.18). 
A more explicit form of the energy density w expressed directly in terms of the field n 
reads 

w=^[n^ + {d^u ■ d^u)] + 1(1 + h^)nl (5.20) 

where we find no trace of the nonrelativistic term h • (n x h) of Eq. (5.16). The energy 
W is now measured in units of 2\ples^ J . 

Therefore the extended nonlinear a model (5.12) may be cast in the standard 
Hamiltonian form 

^a^^. ^a^-YT'^ a = 1,2, 5.21 

OT^a OWa 

18 



which will be the basis for our subsequent theoretical discussion. For instance, the 
conserved linear momentum should be given by 



p^ = - j TTad^i/ja dxdy, (5.22) 

and the angular momentum by 

£ = - y TVaSi^^x^d^^jja dxdy. (5.23) 

Using the canonical variables of Eq. (5.18) a more explicit form of linear momentum 
reads 

P^ = - J{Qd^Q + sin^ G ^d^^ + h sin^ 6 d^^)dxdy (5.24) 

and is in agreement with the expression quoted in Ref. [16]. However the field dependent 
term in Eq. (5.24) leads to an improper integral for vortex configurations for which O ~ ^ 
and $ ~ at spatial infinity. One would think that this difficulty may be resolved by 
modifying the momentum density by a total divergence, thus replacing sin^ 9^$ by 
(sin2e-l)a^$ = — cos G 9^$, which would indeed lead to proper behavior at infinity 
where cos G = 0. Nevertheless the ambiguity would then be shifted to the origin of 
the vortex where is singular, because $ is multivalued, while G = or tt and 
cosG = ±1 0. 

The ambiguities in the linear momentum signal an important link between the 
dynamics and the underlying topological complexity, in analogy with the situation 
previously analyzed for FM bubbles [5, 6]. In turn, the relativistic dynamics of AFM 
vortices studied in earlier sections should be radically altered by the applied field. A 
complete resolution of the ambiguities is given in subsection C after the ground is 
prepared in subsection B. 



B. Vorticity and the stress tensor 

For any 2D field theory that can be brought to the standard Hamiltonian form 
(5.22) one may define a (scalar) vorticity 

7 = ^fJ.iydf.TTadt.lpa, (5.25) 

where e^,^ is the 2D antisymmetric tensor (^n = = £22? £12 = 1 = —£21)- Terminology 
is borrowed from fiuid dynamics because the quantity 7 shares with ordinary vorticity 
several formal properties. The time derivative of 7 is calculated from the Hamilton 
equations to yield 

7 = -£fj,ud^j,Tty (5.26) 

where the vector density 

dW ^ , dW ^ 
r,. = -rrd„i;a + 7 — o^TTa (5.27) 
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is analogous to the "force density" employed by Thiele [17] in the problem of FM 
bubbles. It is not difficult to see that may be written as a total divergence, 

= d\a^x, (5.28) 

where the tensor (7y\ is the stress tensor 

calculated for a specific energy density w. Eq. (5.26) then reads 

7 = -^^udf,dx(Ju\ (5.30) 

and proves to be fundamental for our purposes [5,6]. 

It should be noted that the preceding discussion makes no distinction between 
ordinary field theories and those endowed with nontrivial topological structure or related 
properties. However a clear distinction emerges when we consider the total vorticity 

r = y" ^dxdy = Ef,^ J d^TTadj.ipadxdy, (5.31) 

which is conserved by virtue of Eq. (5.30) for any field configuration with reasonable 
behavior at infinity. One may also write 

^ = £fiu J [d^iTTaduil^a) - TTad^d^il}a]dxdy (5.32) 

to indicate that a vanishing value of the total vorticity is the rule rather than the 
exception. Indeed, under normal circumstances, the first term in (5.32) is shown to 
vanish by transforming it into a surface integral at infinity and the second term also 
vanishes because e^j^dnd,jil)a = for any differentiable function ipa- Yet the above 
conditions may not be met in a field theory with nontrivial topology, a fact closely related 
to the ambiguities discussed in connection with the linear momentum. In general, the 
canonical definition of conservation laws is rendered ambiguous when the total vorticity 
r is different from zero. 

It is then important to examine more closely the definition of vorticity in the current 
model. Substitution of the canonically conjugate fields of Eq. (5.18) in Eq. (5.25) and 
straightforward algebraic manipulation yield the local vorticity 

7 = £tiud^{ii ■ duTo) + kuj, (5.33) 

where 

^ = -^[£p.sin(2e)a,(2e)9^$]. (5.34) 

The first term in Eq. (5.33) is an uncomplicated total divergence which leads to a 
vanishing contribution in the total vorticity F of Eq. (5.31). Thus we may write 



= h J u) dxdy (5.35) 
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and further note that the density uj of Eq. (5.34) resembles the Pontryagin density q 
of Eq. (3.18) except for an overall factor — ^ and the replacement G 20. The latter 
suggests considering the three-component vector N = {Ni, N2, N^) with 

Ni = 277,377,1 = sin(20) cos$, 

N2 = 2n3772 = sin(2e) sin $, (5.36) 
Ns = 2n3-1 = cos(2e), 

which is also a unit vector field (N^ = !)• The density uj may then be written as 

u = —s^^id^N X d^N) ■ N (5.37) 

and should be compared to the standard Pontryagin density of Eq. (3.17). Furthermore 
the field N satisfies the simple boundary condition 

N (0,0,-1), (5.38) 

|X|— >oo 

thanks to the condition 77.3 satisfied by all relevant field configurations, including the 
vortex configurations of Sec. III. The net conclusion is that uj is actually the Pontryagin 
density for the field N and thus yields an integer-valued total vorticity 

T = 27rhK, K = 0,±l,±2,..., (5.39) 

where the integer k will be referred to as the vortex number. Indeed an explicit 
calculation for a single vortex or antivortex (k = ±1) discussed in Sec. Ill confirms 
Eq. (5.39) for any choice of the polarity u, in contrast to Eq. (3.20) that depends on 
both the vortex number and the polarity. 

To complete this level of description of the current model we quote an explicit 
expression for the stress tensor calculated from Eq. (5.29) using as input the energy 
density (5.20). The final result is 

ajyx = w5j,x - {duTo. ■ dxTo), (5.40) 

where the energy density w may be expressed directly in terms of the field n through 
Eq. (5.20). 

Finally we mention that the discussion of this subsection possesses a straightforward 
3D generalization. For instance, Eq. (5.30) becomes 

ii = -Sijkdjdiakh (5.41) 

where Latin indices assume three distrinct values and eijk is the 3D 

antisymmetric tensor. The stress tensor in Eq. (5.41) is obtained by an obvious 3D 
extension of Eq. (5.40) and the vorticity = (71, 72, 73) is given by 

7j = SijkdjT^adki'a = ^ijkdjiii ■ dkn) + huji, (5.42) 
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where the vector density = (a;i, a;2, (^3) reads 

= —SijkidkN X djN) ■ N (5.43) 

and generalizes the scalar Pontryagin density (5.37). Furthermore the vorticity is 
solenoidal (V ■ = 0) essentially by construction and may thus be derived from a vector 
potential ( = V x a) which can be used to define a degree of knottedness or helicity of 
vortex lines, 




(5.44) 



by analogy with the Hopf index in ferromagnets [6] . Such an index may prove to be an 
important issue in the study of 3D antiferromagnets but will not be discussed further 
in the present paper. 



C. Conservation laws 

We now return to the 2D theory and consider the derivation of unambiguous 
conservation laws. Since the main strategy was already explained in the related context 
of FM bubbles [5,6] our description here will address only the essential points adapted 
to the present model. The appearance of a double derivative in the right-hand side 
of the fundamental relation (5.30) suggests that some of the low moments of the local 
vorticity 7 must be conserved. Indeed the linear momentum p = (pi,p2) is given by 

Pfi = -£i^i^L, Iv = J x^^dxdy, (5.45) 
and the angular momentum £ by 



e 



= \ j P^^ c?a;(iy, (5.46) 



where — x'^+y'^. The list of conservation laws is completed by the total magnetization 
H in the third direction, 

fj, = y"[e • (n X h) - hnl]dxdy, (5-47) 

which can be derived directly from the equation of motion. 

The preceding identifications are made plausible by inserting the general expression 
for the vorticity given by Eq. (5.25) in Eqs. (5.45) and (5.46) and by freely performing 
partial integrations to recover the canonical forms of linear and angular momentum 
quoted in Eqs. (5.22) and (5.23) which are plagued by the ambiguities discussed in 
connection with Eq. (5.24). However no such ambiguities occur in Eqs. (5.45) and 
(5.46) because the local vorticity 7 can be obtained directly from the field n, rather 
than the angular variables, and is a particularly well defined quantity; see Eqs. (5.35) 
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and (5.37). In other words, partial integrations should be performed with great care 
and are often unjustified. 

The main point of this theoretical exercise is that the very structure of the 
conservation laws (5.45)- (5.47) suggests a radical change in the dynamical behavior 
of vortices in an applied field {h ^ 0). The effect of a nonvanishing total vorticity 
(r — I'nhK) becomes apparent by considering the transformation of the moments ly of 
Eq. (5.45) under a translation of coordinates x ^ x-|- c where c = (ci, C2) is a constant 
vector, 

+ (5.48) 

which implies a nontrivial transformation of the linear momentum (5.45) when F 7^ 
0. This is surely an unusual property, because linear momentum should be expected 
to remain unchanged under a constant translation of the origin of coordinates, and 
indicates that the moments provide a measure of position rather than momentum. 
Such a fact is made explicit by considering the guiding center vector R = (i?i, R-i) with 
coordinates 

R^ = ^ = - \ Xy^ dxdy, (5.49) 

which transforms as R — > R + c under a constant translation and is thus a measure of 
position of a spin configuration with F 7^ 0. Nevertheless the vector R is conserved. 

A related fact is that the familiar Poisson bracket algebra is significantly affected 
when r ^ 0. Using the canonical Poisson brackets, 

{7r«(x), M^')} = SabSi^ - x'), (5.50) 

and the general expression of the local vorticity (5.25) in the definition of the linear 
momentum (5.45), it is not difficult to establish the relations 

{pi,P2} = r, {R,,R2} = 1/T, (5.51) 

which are strongly reminiscent of the situation in the case of electron motion in a uniform 
magnetic field, the role of the latter being played here by the total vorticity F. 

Similarly the angular momentum (5.46) actually provides a measure of the vortex 
size, a fact made explicit by considering the mean squared radius defined from 

r^ = ^ J{^-Kfldxdy=^-B', (5.52) 

which is also conserved. Needless to say, the conservation laws (5.45) and (5.46) resume 
their ordinary physical significance at vanishing total vorticity (F = 0). 

The observed transmutation in the physical significance of the conservation laws 
of linear and angular momentum implies a radical change in the dynamical behavior 
of topological solitons. For example, a single AFM vortex or antivortex (k = ±1) in 
a uniform magnetic field carries a nonvanishing total vorticity (F = ±27r/j.) and thus 
cannot be found in a free translational motion (R = 0). It is always spontaneously 
pinned or frozen within the antiferromagnetic medium, in contrast to the freely moving 
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vortices occuring in the relativistic theory at vanishing field. Vortex motion can occur 
in the presence of other vortices, but the dynamical pattern is also expected to be 
substantially different from the one obtained at vanishing field in Sec. IV. Specifically, 
interacting AFM vortices should now behave as ordinary vortices in a fiuid or as electric 
charges in a uniform magnetic field, as demonstrated by direct simulations in subsection 
D. 

The preceding descussion was kept deliberately general in order to emphasize that 
the emerging qualitative picture is valid in any field theory for which the total vorticity 
r may be different from zero. However it is now useful to express the conservation 
laws in a more explicit form that takes into account the specific structure of the current 
model. We thus insert the local vorticity of Eq. (5.33) in Eq. (5.45) to obtain the linear 
momentum 

Pi^ = - J[{h - d^n) + hei^^x^(jj]dxdy, (5.53) 

where we have performed a partial integration in the first term which is free of all 
ambiguities. Similarly the angular momentum (5.46) reads 

^ = J [-^i^jy^ni^ ■ 9„n) + ^hp'^uj]dxdy. (5.54) 

As mentioned already, the above conservation laws possess their usual physical 
significance only at vanishing total vorticity, F = 2nhK = 0, which may be achieved 
when either the applied field h or the vortex number k vanishes. Otherwise one must 
consider the guiding center coordinates (5.49) or 



1 



^ 2TThK 



J [-e^j,(h • d„n) + hx^uj\dxdy (5.55) 



and the radius r calculated from Eq. (5.52). It should be noted that the preceding 
conservation laws display some formal similarities to those derived in a model for a 
superconductor [18,19]. 

Finally we return briefly to the 3D theory discussed in the concluding paragraph 
of subsection B and quote the corresponding conservation laws of linear and angular 
momentum 

p = -1 y"(r X )dV, 1 = y"[r X (r X )]dV, (5.56) 

where r = {x,y,z), dV = dxdydz and = (71,72,73) is the vector vorticity field of 
Eq. (5.42). It is interesting that (5.56) are formally idendical to the conservation laws 
derived in fluid dynamics, at least for incompressible fluids; see Eqs. (7.2.5) and (7.2.6) 
ofRef. [20]. 



D. Interacting vortices 

In the presence of a bias field the static vortices of Sec. Ill adjust to a slightly 
different shape. It is not difficult to see that the functional form of a static vortex 
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remains the same as in Eq. (3.13) except that 6 = 9{p) now satisfies the ordinary- 
differential equation 

which differs from Eq. (3.10) only by an additional easy-plane anisotropy with strength 
equal to h^. Consequently Eq. (5.57) reduces to Eq. (3.10) by the simple rescaling 
p = Vl + h^p and the vortex profile 6 = 9{'p) is again given by Fig. 3 with the 
replacement p ^p. More importantly, the auxiliary fields (5.11) or (5.13) now contain 
field dependent terms that are crucial for a correct calculation of the actual spin values 
on the lattice of Fig. 1 or Fig. 2, respectively. 

A pair of like vortices initially at rest is described by the product ansatz (4.2) taking 
into account the field dependent modifications discussed in the preceding paragraph. 
The ensuing time evolution of the vortex pair was obtained numerically. Instead of 
drifting away the two vortices actually begin to rotate around each other, in sharp 
contrast to the situation described in Sec. IV at vanishing field. In other words, each 
vortex moves in a direction perpendicular to the applied force, in analogy with the 
skew defiection of FM bubbles in a field gradient [5,6]. Fig. 9 illustrates the initial 
configuration together with two characteristic snapshots taken at time intervals such 
that the pair had rotated roughly by 90° and 180°, respectively. It is interesting to note 
that a 90° rotation of the pair in real space is always followed by a 90° internal phase shift 
of each vortex. Fig. 10 depicts the actual trajectories obtained by tracking the points 
where |n3| = 1. Inspite of an apparent initial tendency to drift away, the two vortices 
eventually orbit around each other, in complete analogy with the 2D motion of two like 
vortices in an ordinary fluid or two interacting electrons in a uniform magnetic fleld. 
The observed departures of the trajectories of Fig. 10 from a circular shape correspond 
to the well-known Larmor oscillations in the electron problem. These oscillations are 
expected to be smoothed out in the limit of large relative distance. 

The above results are consistent with the qualitative picture suggested by the 
conservation laws of subsection C. The two-vortex system carries a total vortex number 
K = 2 and thus a nonvanishing total vorticity F = Anh. The guiding center calculated 
from Eq. (5.55) is initially located at the origin of the coordinate system and remains 
fixed at all later times. The angular momentum was calculated numerically based on 
Eq. (5.54) and its time evolution is demonstrated in Fig. 11. Although the two pieces 
of Eq. (5.54) acquire a nontrivial time dependence, their sum is fairly well conserved. 
Furthermore the same general picture was obtained by repeating the calculation for 
a pair of vortices with the same vortex numbers (ki = K2) but opposite polarities 
{vi = —1^2)- This result is consistent with the fact that the driving issue is the total 
vorticity F = 27r/i(Ki + K2) which is independent of the polarities, in contrast to the 
ordinary winding number Q of Eq. (3.20). 

The calculation was further repeated for a vortex-antivortex pair initially described 
by the product ansatz (4.8) incorporating the appropriate field dependent modifications. 
Recall that at vanishing field the vortex and the antivortcx arc attracted toward 
each other and are eventually annihilated. The situation is drastically different at 
nonvanishing field. The pair undergoes Kelvin motion roughly along parallel lines that 
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are perpendicular to the line connecting the vortex and the antivortex. Fig. 12 depicts 
the initial configuration together with a snapshot taken at a later instance when the pair 
had moved in formation along the y-ax.is to a distance approximately equal to the initial 
relative separation. Fig. 13 demonstrates the actual trajectories (solid lines) obtained 
by tracking the points where jnaj = 1. Therefore the derived picture is qualitatively 
identical to the Kelvin motion of a vortex-antivortex pair in an ordinary fluid or the 
Hall motion of an interacting electron-positron pair in a uniform magnetic field. 

Returning to the conservation laws we note that the total vorticity of a vortex- 
antivortex pair vanishes {k = hii + hi2 — and hence F = 0). Therefore it is now 
meaningful to interpret (5.53) as the conserved total linear momentum of the system. 
A nonvanishing component develops only along the y-axis, i.e., along the direction of 
motion of the pair, and its conservation is demonstrated in Fig. 14. Again each of the 
two pieces in Eq. (5.53) exhibits a nontrivial time dependence but their sum is fairly well 
conserved. On the other hand, it is still meaningful to define individual guiding centers 
when the pair is widely separated. For example, approximate guiding centers for the 
vortex and the antivortex may be defined by restricting the integration in Eq. (5.55) to 
the right and left half plane and setting k = 1 and k = —1, respectively. The trajectories 
obtained by tracking the above approximate guiding centers are also shown in Fig. 13 
(dashed lines) and are close to two parallel straight lines. The analogy with the motion 
of an electron-positron pair in a uniform magnetic field in now made more definite. The 
actual trajectories of the electron and positron undergo Larmor oscillations along the 
parallel trajectories of their guiding centers, which become increasingly narrower with 
increasing relative separation. The absence of more than one such oscillation in Fig. 
13 is due to our (numerical) inability to follow the motion to a larger distance; see, 
however, a related calculation of interacting vortices in a charged fluid [19]. Finally we 
have verifled that a vortex-antivortex pair (ki + K2 = 0) exhibits Kelvin motion for any 
choice of relative polarities {ui = V2 or vi = —1^2), a fact that reenforces the prominence 
of the total vorticity F in the study of dynamics. 

To summarize, a simple comparison of the results of this subsection to those of Sec. 
IV establishes that the dynamics of AFM vortices is profoundly altered by the applied 
fleld, in remarkable analogy with the familiar Hall effect. Nevertheless one should keep 
in mind that the effect of a bias field on AFM vortices would not have been as drastic 
without the aid of the underlying nontrivial topological structure. 

E. The isotropic antiferromagnet 

The special case of an isotropic Heisenberg antiferromagnet is important for both 
practical and theoretical purposes. The isotropic limit was briefly discussed in the 
concluding paragraphs of Sec. Ill in the absecne of a bias field. It was then mentioned 
that the model possesses mctastable AFM bubbles, instead of vortices, which are 
characterized by the standard Pontryagin index (3.17). However, when a bias field 
is turned on, the picture changes drastically for two reasons. First, the applied field 
itself supplies an effective easy plane anisotropy that leads to vortices instead of bubbles. 
Second, the dynamics of vortices departs significantly from the relativistic dynamics of 
the pure antiferromagnet. 
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In the remainder of this section we shall briefly describe the necessary modifications 
of the formalism to accomodate the isotropic model in a uniform magnetic field. Since 
a single ion anisotropy is no longer available to provide the small parameter e of Eq. 
(2.5), such a parameter is now furnished by the applied field which is assumed to be 
weak: 

e = « 1. (5.58) 

The relevant dynamical equations are then obtained from our earlier results by the 
formal substitution h — > e = (0, 0, 1), or /i — > 1, and by omitting the contribution from 
the single ion anisotropy. Thus the dynamics of the field n is now governed by the 
parameter free Lagrangian 

L = ^[h^ - (a^n • d^n)] + e • (n X h) - ^(e • n)^ (5.59) 

which leads to the equation of motion 

nxf=0, f=h-An + 2(exh) + (n-e)e. (5.60) 
The associated auxiliary fields are accordingly given by 

£ 

m = — p[— (n^ + n^) + (n x h) — n x (n x e)], (5.61) 
2v 2 



for the lattice of Fig. 1, or 



m= ^[(nxh)-nx (n x e)], 

e e 

k=— -Ila;, 1=— -Ily, 



(5.62) 



for the lattice of Fig. 2. Finally we must set h = 1 throughout our discussion of 
conservation laws. 

Therefore the corresponding physical picture can be readily inferred without further 
calculation. Static vortices are formally identical to those of Sec. Ill but their dynamics 
is similar to the nonrelativistic dynamics of the current Sec. V, for any finite value of the 
applied field. In other words, to the extent that topological solitons are relevant for the 
physics of an isotropic antiferromagnet, the dynamical picture is changed significantly 
even by a very weak bias field. 



VI. CONCLUDING REMARKS 



The emphasis in the main text was placed on elementary processes involving only 
two AFM vortices, in order to clearly illustrate an important link between topology 
and dynamics. However further progress in that direction hinges upon the actual 
production of isolated vortices. Although there exist several examples of realistic 
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antiferromagnets that are effectively two-dimesnional, including the parent compounds 
of high-Tc superconductors, there seems to be no direct experimental evidence for 
isolated AFM vortices or bubbles. This situation is in marked contrast to the observed 
abundance of ferromagnetic bubbles, vortices in superfluid helium, or Abrikosov vortices 
in superconductors, and may change in the future. 

The qualitative picture derived from elementary processes must also influence the 
thermodynamics of 2D antiferromagnets. Work in that direction was already presented 
in Refs. [21, 22] for both Ising-like and single-ion anisotropy, while the effect of an 
applied field was considered in Ref. [16]. It is clear that much remains to be done in 
connection with the anticipated Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
which relies on the dynamics of a gas of interacting vortices and antivortices. Suffice 
it to say that the dynamics of vortex-antivortex pairs studied in Sec. IV is radically 
modified by an applied magnetic field discussed in Sec. V. The BKT theory may have 
to be reformulated in a way that clearly refiects the fundamental change of behavior in 
the elementary vortex processes when a field is turned on. 

Perhaps the clearest manifestation of the effect of an applied field will emerge in 
the thermodynamics of an isotropic antiferromagnet. Topological solitons at vanishing 
field are metastable AFM bubbles that obey relativistic dynamics. However the 
smallest external field will trigger Hall dynamics for AFM vortices which become the 
relevant topological excitations. Again a successful BKT theory must reflect this abrupt 
transition in the limit of vanishing fleld. 

Our discussion is concluded with some comments on a variation of the main theme 
that has not been treated in this paper. An easy-axis anisotropy (g < 0) would lead 
to a ground (Neel) state that is polarized along the easy axis. Therefore topological 
solitons must then satisfy the simple boundary condition (3.16) and would be AFM 
bubbles classifled by the standard Pontryagin index (3.17). However an application of 
the Derrick theorem either in its original form or its extended version discussed in the 
Appendix leads to the conclusion that such solitons do not exist either with finite or 
infinite energy. This is a notable difference from FM bubbles that occur in easy-axis 
ferromagnets. The latter are stabilized by a combination of the effects of the long-range 
magnetostatic field created in a ferromagnetic film and of an applied bias field [6] . 

Nevertheless, when an easy-axis antiferromagnet is immersed in a uniform magnetic 
field pointing along the symmetry axis, an effective easy-plane anisotropy is produced 
that competes with the easy-axis anisotropy. And, when the field exceeds a certain 
critical value, a spin-fiop transition takes place from the Neel state polarized in the 
third direction to a canted state of the type shown in Fig. 8 which exhibits azimuthal 
degeneracy. Consequently AFM vortices reappear above the critical value of the applied 
field and their dynamical properties are very similar to those discussed in Sec. V. 
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APPENDIX: VIRIAL THEOREMS 



We derive a slight generalization of the Derrick theorem [7] that is not based on 
the assumption of finite soliton energy. The method is an elementary extension of our 
earlier work on magnetic bubbles [6]. We consider only static solutions for which the 
stress tensor (5.40) in the absence of an external field reduces to 



<7ux = w5,x - [{d^QdxQ) + sin^ e (a,$aA$)], 



w = l[{df,ed^e) + sin^ G + cos^ &], 



(A.l) 



and satisfies the continuity equation 

dxCT.x = 0, (A.2) 

on account of the static Hamilton equations (3.5). 

A series of virial relations may now be derived by taking suitable moments of (A.2), 
the simplest possibility being 

/ x^dxafj^xdxdy = 0, (A.3) 
J A 

where the integration extends over some finite region A in the 2D plane. An application 
of the divergence theorem yields 

/ (Tfj,ydxdy = i x^a/j^xdsx, (A.4) 

J A Jc 

where the integral in the right-hand side is taken over the closed curve C surrounding the 
area A; by convention the line-element vector {dsi,ds2) is perpendicular to the curve. 
On the assumption that all fields approach asymptotic values such that the volume 
integrals in the left-hand side are finite and the surface integrals in the right-hand side 
vanish, when the integration extends to infinity, one obtains the virial relations 

a^^dxdy = 0, (A. 5) 

where n^u — 1 or 2 are taken in any combination. It is not difficult to see that relations 
(A. 5) may also be obtained by a Derrick-like argument where coordinates are rescaled 
according to x ^ A11X + A12?/, y — > )^2iX + X22y and we further demand that the resulting 
total energy be stationary at An — 1 — A22 and A12 = = A21. In fact, the original 
scaling relation of Derrick may be recovered by taking the trace of Eq. (A. 5): 



/ 



tradxdy = 0. (A. 6) 



Substitution in Eq. (A. 6) of the explicit stress tensor (A.l) reproduces Eq. (3.7) which 
contradicts the existence of solitons with finite energy, the latter requirement being 
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implicit in the transition from Eq. (A. 4) to Eq. (A. 5) as well as in the original argument 
of Derrick [7] . 

In order to probe the existence of interesting static solutions with infinite energy, 
we return to Eq. (A. 4) and take the limit A ^ oo more carefully. For definiteness let us 
restrict our attention to the axially symmetric configuration of Eq. (3.9), i.e., = 6{p) 
and $ = K{(j) + (l)o) with n = ±1, which is inserted in the stress tensor (A.l) to yield the 
reduced elements 

(711 = -Ucos{2(p) + V, (722 = Ucos{2(p) + V, 

(A.7) 

(712 = -U sin(20) = (721, 



where the functions 



U = 



dp 



V=]- cos^ 9, 
2 



(A. 



depend only on the radial coordinate p. Furthermore the area A in Eq. (A. 4) is taken to 
be a circle of radius R and the Cartesian components of the line-element vector ds are 
expressed in terms of the polar coordinates as c^si — Rcoscpdcp and ds2 = RsiiKp d(l). 
Relations (A. 4) applied for z/ = 1 or 2 taken in all combinations reduce to the single 
virial relation 

R 

27r pdpV{p) = 7rR^[V{R)-U{R)], (A.9) 

which must be satisfied for any radius R. Indeed, if we use the numerically calculated 
vortex profile 9 = 9{p) shown in Fig. 3 to calculate the functions U = U{p) and V = V{p) 
from Eqs. (A. 8), relation (A.9) is verified for all R. In particular, the limit i? — > oo yields 
a nonvanishing contribution in the right-hand side thanks to the centrifugal term in the 
function U which leads to U{R) ~ —1/2R^ because 6* — > ^. The net result of this limit 
is 

I 2n pdp cos^ 9 = 7r, (A. 10) 

Jo 

which establishes the virial relation announced in Eq. (3.8). This relation no longer 
contradicts the existence of nontrivial vortex solutions but actually predicts analytically 
their anisotropy energy. However the exchange energy is (logarithmically) divergent 
thanks again to the centrifugal term. 
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FIGURE CAPTIONS 



Figure 1 : Illustration of the dimerization process for a finite portion of the square lattice 
cut along the diagonals. Solid and open circles denote the sites of the two intertwining 
sublattices. 

Figure 2: Illustration of the tetramerization process of a finite portion of the square 
lattice. 

Figure 3: The numerically calculated profile of a static vortex at weak anisotropy. 

Figure 4: Spin vectors of a vortex and an antivortex projected on the (12) plane. 
Spins are shown at every fourth site of the original lattice and distances on the grid are 
measured according to Eq. (2.24) applied with £ = 0.1. 

Figure 5: Three characteristic snapshots of a head-on collision of two like vortices 
originating at a relative distance d = 5 on the a;-axis with initial velocities vi = (0.65, 0) 
and V2 = (—0.65, 0). After collision the two vortices scatter at 90° and suffer an internal 
phase shift also equal to 90°. Vectors represent the projection of the field n on the (12) 
plane. 

Figure 6: Three characteristic snapshots of the annihilation process of a vortex and 
an antivortex that are initially at rest at a relative distance d = 5 on the x-axis. The 
vortex and the antivortex converge toward each other and are eventually annihilated 
into spinwaves. Vectors represent the projection of the field n on the (12) plane. 

Figure 7: Level contours of the energy density corresponding to the vortex-antivortex 
annihilation process shown in Fig. 6. After collision two distinct energy lumps emerge 
along the positive and negative j/-axis but eventually dissipate into spinwaves. 

Figure 8: Schematic illustration of the ground state in the presence of an external 
uniform magnetic field along the symmetry (third) axis. The canting angle 5 is given 
by Eq. (5.3). 

Figure 9: Time evolution of a pair of like vortices in the presence of a bias field 
h = 1. The two vortices are initially at rest, at a relative distance d = 5 on the x-axis, 
and subsequently orbit around each other instead of drifting away. Together with the 
initial configuration of the field n projected on the (12) plane (upper entry) we provide 
snapshots at instances when the pair had rotated by 90° (middle entry) and 180° (lower 
entry). A 90° rotation in real space is always followed by a 90° internal phase shift. 

Figure 10: Trajectories of the two rotating vortices of Fig. 9 obtained by tracking the 
points where |n3| = 1. The two vortices originated at points A and B on the x-axis and 
the process was interrupted after the pair had completed a 180° rotation. 

Figure 11: Evolution of the angular momentum of the vortex pair of Figs. 9 and 
10 calculated from Eq. (5.54). The total angular momentum £ (solid line) is fairly 
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well conserved, except at late times when numerical instabilities develop on the finite 
lattice. The lower (upper) dashed line depicts the nontrivial time dependence of the 
first (second) term in Eq. (5.54). 

Figure 12: Evolution of a vortex-antivortex pair in the presence of a bias field 
h = 1. The vortex and the antivortex are initially at rest at the points (2.4, —2.4) 
and (—2.4, —2.4) of the xy-plane (upper entry). Instead of converging toward each 
other and annihilating, the pair moves in formation along the y-axis (Kelvin motion) as 
demonstrated by the snapshot of the lower entry taken at a later instance (r ~ 20). 

Figure 13: Trajectories of the vortex and the antivortex of Fig. 12 originating at the 
points B = (2.4, —2.4) and A = (—2.4, —2.4). The actual trajectories (solid lines) were 
obtained by tracking the points where |?i3| = 1. The trajectories of the guiding centers 
(dashed lines) were calculated from Eq. (5.55) restricted to the corresponding half 
planes and are remarkably close to two parallel straight lines. 

Figure 14: Evolution of the linear momentum of the vortex-antivortex pair of Figs. 
12 and 13 calculated from Eq. (5.53). (We display only the y-component because the 
x-component vanishes.) The total linear momentum (solid line) is fairly well conserved, 
whereas the lower (upper) dashed line depicts the nontrivial time dependence of the 
first (second) term of Eq. (5.53). 
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